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A  zero-dimensional  time-dependent  high-pressure  slab  capillary  discharge  model  with 
capillary  wall  thermal  conduction  and  wall  radiation  absorption  is  presented.  The  model 
includes  a  resistor-inductor-capacitor  circuit  and  a  heat-transfer  radiation  model  based  on 
radiation  database  constructed  using  PrismSPEC,  a  commercially  available  radiation 
software,  to  calculate  the  radiation  heat  flux  output  from  a  uniform  plasma  slab.  The  model 
also  includes  a  model  of  the  thin  transition  boundary  layer  between  the  uniform  plasma  core 
and  the  ablative  capillary  walls.  This  transition  boundary  layer  model  is  used  to  obtain  the 
boundary  conditions  connecting  the  plasma  core  parameters  with  the  parameters  at  the 
ablative  surface  to  calculate  the  thermal  and  radiation  heat  fluxes  at  the  capillary  walls.  The 
radiation  wall  absorption  coefficient  and  the  wall  thermal  conductivity  are  assumed  to  be 
constants  independent  of  wavelength  and  wall  temperature.  Thus,  the  model  self-consistently 
calculates  plasma  parameters  and  distribution  of  wall  temperature  vs.  time.  The  model  can 
be  used  to  model  high-pressure  ablative  capillary  discharge  for  plasma  thrusters. 


Nomenclature 

A  =  average  mass  of  heavy  particles  (neutrals  and  ions)  in  a  atomic  units 

C  =  capacitor  of  electrical  circuit  [F/m] 

Cs  =  sound  speed 

Ca  =  carbon  fraction  of  polyethylene  molecule 

cP  =  specific  heat  at  constant  pressure 

Da  =  slab  capillary  gap 

d  =  width  of  the  initial  wall  temperature  distribution 

end  =  index  denoting  condition  at  the  exit  plane 

Fmd  =  radiation  heat  flux  incident  on  the  transition-boundary  layer 

Ha  =  hydrogen  fraction  of  polyethylene  molecule 

//„  =  thickness  of  the  capillary  wall 

hCH  =  enthalpy  of  a  CH-molecule  to  bring  it  from  the  capillary  wall  to  the  plasma  volume 
I  =  discharge  current  per  unit  of  slab  length 

kB  =  Boltzmann  constant 

L  =  inductor  of  electrical  circuit  [H  m] 

La  =  capillary  length 

Mc  =  mass  of  a  hydrogen  atom  in  atomic  units 

Mh  =  mass  of  a  carbon  atom  mass  in  atomic  units 

M0  =  unit  of  atomic  mass 

n  =  normalized  gas  number  density,  a  variable  of  the  Knudsen  layer  model 

na  =  number  density  of  neutral  atoms 

nCH  =  molecular  number  density  of  ablated  molecules 

n,  =  ion  number  density 

P  =  pressure 
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Q  =  charge  of  capacitor  [Q/m] 

R  =  resistor  of  electrical  circuit  [Q-m] 

Rpiasma  =  plasma  resistance  [Q-m] 

sat  =  index  denoting  vapor  saturation  conditions 

sur  =  index  denoting  conditions  for  ablative  wall  surface 

T  =  plasma  temperature 

Th  =  maximal  wall  surface  temperature 

Tsur  =  temperature  of  the  ablative  surface 

Tw  =  wall  temperature 

t  =  time 

u  =  normalized  directed  velocity,  a  variable  of  the  Knudsen  layer  model 

Vt  =  thermal  gas  velocity 

V  =  normalized  thermal  velocity,  a  variable  of  the  Knudsen  layer  model 

x  =  axis  directed  across  the  slab 

Z  =  average  ion  charge  state 

p  =  ratio  of  thermal  pressure  to  magnetic  pressure 

P  =  variable  of  the  Knudsen  layer  model 

Y'ai  =  molecular  ablation/absorption  rate 

r cH.max  =  maximal  molecular  ablation  rate 

y  =  ratio  of  specific  heats 

A  So  =  sum  of  vaporization,  dissociation,  ionization  and  electronic  excitation  energies 

Sch  =  inner  plasma  energy  per  CH  molecule 

?7  =  radiation  absorption  coefficient  of  the  wall  material 

©  =  energy  flux  incoming  into  the  capillary  wall 

kw  =  wall  material  specific  heat 

fl  =  radiation  heat  source  in  the  bulk  of  capillary  wall 

o  =  electrical  conductivity 

pw  =  wall  material  mass  density 

<p  =  ionization  ratio 

XK  =  wall  thermal  conductivity 


I.  Introduction 

The  role  of  space-based  systems  for  both  commercial  and  government  customers  continues  to  evolve, 
continually  providing  new  requirements  for  the  development  of  satellite  propulsion  systems.  As  it  is  very  efficient 
(50-70%  kinetic -to-electric-power  ratio),  high-pressure  (5  106-5  108  Pa)  ablative  capillary  discharge  is  a  good 
candidate  for  high-power  (tens  kilowatts)  plasma  thrusters. 

A  schematic  of  a  slab  capillary  discharge  thruster  is  shown  in  Fig.  1.  The  discharge  maintains  a  resistive  arc 
through  a  narrow  insulating  capillary  by  the  continual  ablation  of  an  injected  mass  or,  more  commonly,  the  capillary 
wall  material  as  shown  in  Fig.  1.  The  ablative  capillary  discharge  can  be  described  via  a  three-layer  configuration. 
The  outermost  layer  is  a  solid  wall,  usually  some  form  of  polyethylene,  occasionally  Teflon®  or  some  other 
insulating  material.  Material  from  the  wall  is  evaporated,  entering  the  thin  transition  boundary  layer,  where  it  is 
dissociated,  ionized  and  heated  to  plasma  temperature.  The  innermost  layer  is  the  plasma  core.  The  closed  end  of  the 
capillary  (left  side  of  Fig.  1)  is  one  electrode;  the  other  electrode  (right  side  of  Fig.  1)  is  an  open  end  through  which 
the  plasma  can  flow  and  expand.  The  polarity  has  no  appreciable  effect  on  the  discharge  operation.  Ohmic  heating  is 
responsible  for  plasma  heating,  ionization,  ablation,  and  radiation.  In  the  model,  the  sonic  condition  (Ma  =  1  in  Fig. 
1)  is  assumed  to  exist  at  the  open  end  of  the  capillary. 

In  terms  of  plasma  propulsion  concepts,  the  capillary  plasma  thruster  must  satisfy  two  main  conditions.  First, 
the  capillary  plasma  should  be  in  local  thermodynamic  equilibrium  (LTE),  indicating  that  all  sorts  of  particles, 
electrons,  ions  and  neutrals,  have  the  same  temperature.  Then,  the  electrical  current  will  heat  both  the  electrons  and 
the  heavy  particles  (i.e.,  not  only  the  electrons),  thus  providing  high  thrust  per  watt.  Second,  as  it  is  well  known, 
pinching  the  capillary  discharge  leads  to  a  formation  of  a  narrow  plasma  core  region  with  a  high  electron 
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temperature  and  a  large,  relatively  cold  peripheral  plasma  region.  This  also  decreases  the  efficiency  of  the  capillary 
thruster  because  only  a  small  amount  of  input  discharge  energy  is  transferred  to  the  heavy  particles  and  almost  all  of 
the  input  energy  is  transferred  to  the  electrons.  To  minimize  the  pinching  of  the  capillary  discharge,  the  ratio  of  the 
thermal  pressure  to  the  magnetic  pressure,  /?,  in  the  capillary  plasma  thruster  should  be  greater  than  unity. 

The  first  paper  to  discuss  the  use  of  ablative  capillary  discharge  for  propulsion  was  published  in  the  United 
States  in  19821.  In  it,  the  authors  introduced  their  concept  of  the  pulsed  electrothermal  thruster;  this  group  later  built 
and  tested  such  devices2' 3.  Since  then  many  theoretical  papers  have  been  published  exploring  the  applications  of 
such  discharges  to  electrothermal  capillary  thrusters,  electrothermal  launchers,  and  electrothermal  chemical  guns;  a 
review  of  these  papers  with  the  corresponding  references  can  be  found  in  a  recent  theoretical  paper4.  It  should  be 
noted  that  in  this  paper4  the  author  for  the  first  time  used  a  radiation  heat-transfer  model  based  on  a  radiation 
database.  This  database  was  constructed  using  PrismSPECT,  a  commercially  available  radiation  software5  that 
calculates  the  radiation  spectrum  output  from  a  uniform  plasma  slab.  Thus,  unlike  earlier  models,  this  model  does 
not  use  any  asymptotic  radiation  models,  such  as  blackbody  radiation  models  with  a  grey  factor  or  volumetric 
bremsstrahlung  (free-free)  radiation  model,  but  self-consistently  calculates  the  radiation  heat  flux  at  the  transition 
boundary  layer.  Fig.  1 .  However,  it  should  be  stressed  that  model4  neglects  the  radiation  heating  of  the  bulk  capillary 
wall,  assuming  that  all  radiation  entering  the  transition  boundary  layer  is  absorbed  into  this  region  and  expended  on 
the  ablation  of  wall  material.  On  the  other  hand,  as  it  has  been  indicated  in  works6' 7' 8  the  heat  losses  in  the  capillary 
wall,  can  play  a  noticeable  role  in  heat  transfer  of  relatively  high-pressure  (2T05  -  2 TO6  Pa)  capillary  discharges. 
Therefore,  extending  model4  to  the  case  of  heat  transfer  in  the  capillary  wall  is  important  to  better  understand  the 
physics  of  high-pressure  capillary  discharge  and  engineering  aspects  of  the  capillary  discharge  thrusters. 

The  present  paper  extends  the  model4  to  the  case  of  heat  transfer  in  the  bulk  of  capillary  wall  by  taking  into 
account  the  capillary  wall  thermal  conduction  and  radiation  absorption.  A  description  of  the  model  and  numerical 
results  are  presented  in  Sections  II  and  III,  respectively,  and  conclusions  are  given  in  Section  IV. 


II.  Description  of  the  Model 


A.  Assumptions 

The  following  assumptions  are  made  in  the  model: 

1.  The  temperatures  of  electrons  and  heavy  particles  (ions  and  neutrals)  are  equal,  that  is,  the  LTE 
condition  is  achieved. 

2.  Magnetic  pressure  is  much  smaller  than  thermal  pressure,  that  is,  J3  »1 . 

3.  The  conduction  heat  flux  in  the  plasma  core  region  is  much  smaller  than  the  radiation  flux  at  the  inner 
boundary  of  the  transition  boundary  layer  and  can  be  ignored. 

4.  Plasma  composition  can  be  calculated  using  the  Saha  equation,  that  is,  local  plasma-chemistry 
equilibrium  is  achieved. 

5.  Heating  of  the  plasma  by  the  viscosity  drag  force  is  small  and  can  be  omitted. 

6.  The  capillary  wall  is  polyethylene,  C4H9. 

7.  The  capillary  is  a  4  mm  slab. 

8.  The  gas  in  the  plasma  core  region  and  in  the  boundary  transition  layer  is  fully  dissociated. 

9.  The  temperature  of  the  ablative  gas  entering  the  plasma  core  region  from  the  transition  boundary  layer. 
Fig.  1,  is  equal  to  T ,  the  average  temperature  of  the  capillary  plasma. 

10.  The  boundary  transition  layer  is  assumed  to  be  thin  and  is  treated  in  the  model  as  the  interface  between 
the  plasma  core  and  the  capillary  wall. 

All  ten  of  these  assumptions  have  been  analyzed  in  detail  in4  and  will  not  be  discussed  in  this  paper. 
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11.  Unlike  in4,  where  the  heat  losses  into  the  capillary  wall  were  neglected,  in  this  paper  we  take  them  into 
account.  We  assume  that  the  capillary  wall  radiation  absorption  coefficient,  77,  is  a  constant, 
independent  of  wavelength  and  wall  temperature.  This  assumption  is  a  stretch  because  the  extinction 
coefficient  of  polyethylene  is  considerably  dependent  on  wavelength9'  10  and  temperature.  However, 
since  our  model  is  not  capable  of  calculating  the  radiation  spectrum  at  the  capillary  wall  anyway,  and, 
to  the  best  of  our  knowledge  currently  there  is  no  experimental  or  theoretical  data  on  optical  properties 
of  polyethylene  as  a  function  of  temperature,  we  have  used  77  as  an  input  parameter. 

12.  The  thermal  conduction  coefficient  of  polyethylene  is  a  constant,  independent  of  temperature. 

In  our  model  we  use  radiation-plasma-composition  database4.  This  database  was  constructed  using 
PrismSPECT  software5  for  C4H9  capillary  wall  composition  and  described  in  detail  in4. 

B.  Basic  equations 

The  equations  describing  mass  and  energy  conservation  laws  in  the  plasma  core  region  of  a  slab  capillary 
discharge.  Fig.  1,  can  be  written  as 


D„  ■  L„ 


dn. 


CH 


dt 


-2  La-TCH  ~Da-  nCH  end  ■  Cs  end 


(1) 


D„  •  L„ 


d(nCH  ■ SCH ) 
dt 


1 2  •  L„ 
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s.end 


7  'k  B  • Tend  '(l  +  Z) 

y* 

A  — 

r  Ca-Mc+Ha-MH) 

A-M0 

,  /l  — 
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(2) 

(3) 

(4) 

(5) 

(6) 

(7) 


Ca  =  4,  Ha  =  9,  y=  5/3,  cP  =  y/(y-  1),  and  <p=  77/(77 ,  +  na)  is  the  ionization  ratio;  the  vaporization  and  dissociation 
energies  are  considered  as  constants  in  the  model  and  as  input  parameters  in  the  code;  a  detailed  explanation  of  how 
electron  states,  plasma  conductivity,  and  plasma  composition  are  calculated  can  be  found  in4;  the  index  “end” 
denotes  conditions  at  the  exit  plane  and  the  absence  of  this  index  denotes  conditions  averaged  over  the  volume  in  the 
“middle”  of  capillary.  In  the  model  we  use  isentropic  flow  relations 1 1  and  assume  that  the  average  ion  charge  and 
ionization  ratio  are  preserved  through  the  capillary  up  to  the  exit  plane,  Zmd  =  Z,  and  <pend  =  <p  . 

In  the  mass  conservation  equation,  Eq.  (1),  the  first  term  in  the  right-hand  side  describes  the  flow  of  “CH- 
polyethylene  molecules”  arriving  to  the  plasma  core  region  from  the  ablative  surface,  and  the  second  term  is  the 
flow  of  CH-polyethylene  molecules  leaving  the  capillary  chamber  though  the  open  capillary  end.  It  is  worth  noting, 
that  the  case  of  rCH  <  0  corresponds  to  the  plasma  decay  case  where  the  particles  are  absorbed  by  the  capillary  wall 
(not  ablated  from  the  wall);  in  Eq.  (1-2)  the  mass  flux  is  directed  to  the  capillary  chamber.  In  Eq.  (2)  the  first  term  in 
the  right-hand  side  is  Ohmic  heating  of  the  plasma,  the  second  term  is  enthalpy  flow  leaving  the  capillary  with  the 
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plasma  jet  through  the  open  capillary  end,  and  the  third  term  is  energy  losses  in  the  bulk  of  the  capillary  wall.  It 
should  be  stressed  that  we  assume  that  the  heat  flux  can  be  directed  only  from  the  plasma  to  the  wall,  Eq.  (3).  This 
has  a  perfect  physical  explanation:  in  high-pressure  capillary  discharges  the  temperature  of  the  plasma  core  region  is 
always  larger  than  the  temperature  of  the  polyethylene  capillary  which  is  usually  much  less  than  0.1  eV.  As  one  can 
see,  taking  0  =  0,  i.e.,  no  heat  losses  in  the  bulk  of  the  capillary  wall,  we  recover  model4. 

The  heat  transfer  in  the  bulk  of  the  capillary  wall  can  be  written  as 


'  Pw 


°r..-Xw.£Lf..n 


dt 


dx z 


(8) 


where 


f  n  ■  exp (-77  •  X)  •  (Frad  -  r CH  hCH) 


n  = 


1 1  ■  exp(-t7  •  X)  ■  Frad 


0 


if  {Frad  >  FCH  '  hCH  )&(Tcff  >0) 
if  (rCH  <  o) 
f  iFrad  <  FCH  '  hCH  ) 


(9) 


The  first  case  in  Eq.  (9)  corresponds  to  the  case  of  the  wall  ablation  process.  In  this  case  the  radiation  leaving  the 
plasma  core  region.  Fig.  1,  ablates  the  capillary  wall,  dissociates,  ionizes,  and  heats  the  ablative  vapor  up  to  the 
plasma  temperature;  it  is  then  absorbed  in  the  bulk  of  capillary  wall;  and  finally,  it  leaves  the  capillary.  The  second 
case  corresponds  to  the  plasma  decay  process,  where  the  mass  flux  is  directed  from  the  capillary  chamber  to  the 
capillary  wall.  In  this  process  all  the  plasma  radiation  is  absorbed  by  the  bulk  of  capillary  wall  or  leaves  the 
capillary;  the  enthalpy  flux  of  particles,  TCh  hcH ,  is  absorbed  by  the  wall  capillary  surface  and  is  taken  into  account 
in  the  boundary  conditions  at  the  inner  surface  of  the  capillary  wall,  Eq.  (20)  case  2.  The  third  case  in  Eq.  (9) 
corresponds  to  the  situation  were  all  radiation  incoming  into  the  boundary  transition  layer  from  the  plasma  core 
region  is  absorbed  in  it  and  expended  for  ablation  of  the  wall,  the  case  of  no  heat  looses  in  the  bulk  of  the  capillary 
wall. 

It  should  be  stressed  that  in  our  model  we  do  not  take  into  account  the  change  in  the  slab  gap,  Da ,  during  the 
discharge:  the  slab  gap  in  our  radiation  database  is  assumed  to  be  constant  and  cannot  be  varied.  This  limitation  is 
not  critical  for  short  pulse  capillary  discharges;  in  our  calculation  the  increase  in  the  Da  was  about  a  few  microns. 
However,  for  long  pulse  discharges  the  effect  of  changing  the  capillary  radius  due  to  ablation  has  to  be  taken  into 
account. 

In  our  model  we  employed  the  resistor-inductor-capacitor  circuit.  Fig.  2,  to  self-consistently  calculate  the 
parameters  of  the  capillary  discharge  and  the  temperature  distribution  in  the  capillary  walls  vs.  time.  The  equation 
describing  the  circuit  can  be  written  as 
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(10) 


The  initial  parameters  for  this  circuit  are  initial  capacitor  voltage  and  initial  circuit  current  I  =  dQ/dt. 

C.  Boundary  conditions 

Boundary  conditions  connecting  the  plasma  core  parameters  with  the  parameters  at  the  ablative  surface  are  as 
follows: 
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where  Eq.  (11)  calculates  the  flux  of  the  ablative  polyethylene  molecules;  Eq.  (12)  states  that  the  pressure  in  the 
plasma  core  is  equal  to  the  vapor  pressure  in  the  Knudsen  layer  at  the  ablative  surface  (the  Knudsen  layer  is  a  part  of 
the  transition  layer);  and  Eqs.  (13)  -  (15)  describe  the  Knudsen  layer  at  the  ablative  surface  using  Anisimov’s 
velocity  distribution  function12;  here  we  have  adopted  the  Knudsen  layer  model13.  In  the  model  we  use  the  following 
expressions  for  the  equilibrium  polyethylene  vapor  pressure  Psat  1 ,  thermal  velocity  of  the  evaporated  atoms  VTsun 
and  density  nCH.sat, 
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In  Eq.  (18)  we  have  assumed  that  all  vapor  leaving  the  ablative  surface  is  fully  dissociated.  It  worth  noting  that  in 
this  Knudsen  layer  model,  we  ignore  the  thermal  conductivity  in  the  Knudsen  layer,  assuming  that  radiation  heat 
flux  at  the  ablative  surface  is  much  larger  than  thermal  conduction  [14], 

Boundary  conditions  at  the  outer  and  inner  boundaries  of  the  capillary  walls  are 
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where  xmr  =  0  and  x  =  Hw  are  the  coordinates  of  the  inner  and  outer  boundaries  of  the  capillary  walls.  Eq.  (19) 
specifies  no  conduction  heat  flux  at  the  outer  boundary  of  the  capillary  wall.  The  first  case  in  Eq.  (20)  specifies  no 
conduction  heat  flux  at  the  inner  boundary  of  the  capillary  wall  in  the  wall  ablation  regime.  The  second  case 
corresponds  to  the  absorption  wall  regime,  the  regime  in  which  the  mass  flux  is  directed  from  the  plasma  to  the 
capillary  walls.  In  this  regime  the  incident  particles  bring  their  enthalpy  to  the  capillary  wall  surface,  inducing  the 
corresponding  conduction  heat  flux  at  the  capillary  wall.  The  third  case  specifies  the  surface  wall  temperature  as  a 
maximum  wall  temperature  Th  that  can  provide  the  maximum  ablation  rate  due  to  radiation: 
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Thus,  Eq.  (20)  as  well  as  Eqs.  (3)  and.  (9),  state  that  heat  flux  can  be  directed  only  from  the  plasma  to  the  wall  and, 
therefore. 


1  Cl l  -  1  a/.mix 


(22) 


D.  Algorithm 

Let  us  now  describe  the  algorithm  for  calculating  capillary  discharge  parameters  and  the  temperature 
distribution  in  the  capillary  wall: 

1.  Assume  initial  values  of  nah  and  T,  and  r„,(x),  /,  and,  and  Q  . 

2.  Calculate  all  plasma  parameters  for  given  nab  and  T. 

3.  Calculate  V CHmax  busing  Eq.  (21). 

4.  Calculate  the  maximum  allowed  wall  surface  temperature  Tb  by  substituting  TCHmax  and  P (nCH,T)  into 
the  system  of  Eqs.  (1 1)  -  (15). 

5.  Calculate  rCHby  substituting  P (nCH,T)  and  Tw(xsur)  into  the  system  of  Eqs.  (1 1)  -  (15). 

6.  Determine  the  boundary  conditions  at  x  =  xsur,  Eq.  (20): 

6a.  if  ( Tw(xsur )  <Tb)  &  (Fch>  0)  then  apply  case  1; 

6b.  if  ( TyXxsur )  <  Tb)  8l(Tch<  0)  then  apply  case  2; 

6c.  if  TJxsur)  >Th  then  apply  case  3  and  set  T a,  equal  to  T cu.mm', 

7.  Calculate  new  time:  t  =  t+  r,  where  t  is  the  initial  time  and  r  is  the  time  step. 

8.  Calculate  new  nCH,  and  7’  by  solving  Eq.  (1)  and  (2): 
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9.  Calculate  new  TJx)  by  solving  Eq.  (8)  with  boundary  conditions  (20). 

10.  Calculate  new  Q  and  I  solving  Eq.  (10): 

Qt+T  =  Qr  - 11  ■  r 


jt  +  T  _  jt 


1  ~T 


( Rpla  sma  +  R) 


+  r  ■ 


Qr 

CL 
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(26) 


11.  Go  to  step  1. 

In  this  work,  we  have  used:  the  iteration-divide-segment-by-two-method  to  solve  Eq.  (24)  and  the  system  of  Eq. 
(11)-  (15);  a  standard  implicit  scheme  to  solve  Eq.  (8);  and  a  non-uniform  geometric-progression  grid. 
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IV.  Numerical  results 


In  the  numerical  results  presented  in  this  section,  the  combined  energy  of  evaporation  and  dissociation  was 
taken  as  50  eV  per  C4-H9  polyethylene  molecule,  the  capillary  length  was  chosen  as  8  cm,  Da  was  4  mm,  the 
parameters  of  circuit  were  L  =  4.661-10'8  H  m,  R  =  1.4546  Q-  m,  C  =  1.7507  Q/m,  Fig.  1,  the  initial  voltage  on  the 
capacitor  and  the  initial  discharge  current  were  2500  V  and  1.496- 103  A/m.  The  initial  temperature  distribution  of 
the  wall  was  chosen  as, 


rH,(x)  =  (527.56-300.0)- 


Exp[-(x/r/)2] 

1  -  Exp[-(H  w/rj)2] 


^  300.0-  527.56  ■  Exp[-(Hn, /q)2] 
1  -Exp[-(Hw/n)2] 


(27) 


where  a  =  0  corresponds  to  the  ablative  surface,  a  =  //„,  to  the  outer  boundary  of  the  capillary  wall,  and  the  initial 
parameters  of  the  plasma  core  were  chosen  as  nCH  =  1.714- 1023  m"3  and  T  =  1.2555  eV.  To  investigate  how  the 
radiation  absorption  of  the  wall  affects  the  capillary  discharge  we  performed  calculations  for  //  =  0.02,  0.1,  and  0.5 
mm.  It  is  worth  noting,  that  the  selected  initial  conditions  correspond  to  the  parameters  of  the  stationary  capillary 
discharge  with  no  heat  losses  in  the  bulk  of  the  capillary  walls. 

Figs.  3-4  show  the  distributions  of  the  LTE  parameter4  (which  is  the  ratio  of  the  average  exchange  energy 
between  the  electrons  and  the  heavy  particles  due  to  their  collisions  to  the  average  energy  that  an  electron  gains  from 
the  electric  field  between  collisions)  and  /?  vs.  time  for  different  77.  As  one  can  see,  these  parameters  are  much  larger 
then  one;  meaning  that  assumptions  1  and  2  are  well  satisfied  in  our  calculations.  In  the  model  we  have  also  checked 
assumptions  3-5  employing  the  methods4;  they  were  also  well  satisfied  in  our  calculations. 

Now  let  us  analyze  Figs.  5  -  9.  It  is  obvious  that  as  the  radiation  wall  absorption  coefficient  77  increases,  the 
temperature  of  the  ablative  surface  decreases,  as  shown  at  the  positive  slopes  of  Tmr  in  Fig.  5,  leading  to  a  smaller 
vaporization  rate  of  the  wall  material.  This  results  in  smaller  plasma  densities  and  a  correspondingly  smaller  plasma 
pressure  for  larger  77,  as  shown  in  Figs.  6  and  7,  but  in  larger  plasma  temperatures  due  to  the  Ohmic  plasma  heating 
(Fig.  8).  On  the  other  hand,  because  of  the  highly  nonlinear  dependence  of  the  plasma  electrical  resistivity  on  the 
plasma  temperature  and  plasma  composition,  the  discharge  current  does  not  exhibit  a  significant  dependence  on  77, 
as  shown  in  Fig.  9;  however,  it  is  clear  that  with  an  increase  in  77,  the  peak  current  rises.  The  dip  in  nCH  and  the  peak 
in  T  at  the  early  stage  of  the  discharge  seen  in  Figs.  (6)  and  (8)  can  be  explained  by  a  large  imbalance  between  the 
ablation,  rapid  Ohmic  heating,  and  plasma  exhaust  through  the  open  capillary  end  at  the  initial  stage  of  the  capillary 
discharge  development;  the  characteristic  time  of  this  stage  can  be  estimated  as  a  few  LJCS  ,  which  translates  to 
about  20-30  microseconds  in  our  case. 

At  the  negative  slopes  of  Tsur ,  Fig.  5,  all  parameters  of  the  discharge  become  weakly  dependent  of  77,  Figs.  (3-4) 
and  (6-9).  This  observation  can  be  explained  by  the  fact  that  when  Tsur  decreases  with  time  the  heat  losses  in  the 
bulk  of  the  capillary  wall  become  to  be  equal  to  zero  in  the  model.  Fig.  (10).  These  figures  also  demonstrate  that  the 
energy  losses  in  the  capillary  bulk  dramatically  increase  with  an  increase  in  the  radiation  wall  absorption 
coefficients,  as  of  course  we  expect. 

Figs.  11-13  illustrate  wall  temperature  distributions  at  different  moments  of  time  and  for  different  77.  Since  the 
thermal  conductivity  of  the  polyethylene  is  small,  the  characteristic  widths  of  the  temperature  distributions  depend 
on  77.  However,  when  the  ablative  surface  is  cooled  by  the  vaporization  process,  the  length  of  the  positive 
temperature  slopes  depends  on 


V.  Conclusions 

We  have  presented  a  zero-dimensional  model  of  high-pressure  slab  capillary  discharge  model  with  capillary 
wall  thermal  conduction  and  radiation  absorption.  The  model  includes:  a  heat-radiation  model  based  on  radiation 
database  to  self -consistently  calculate  the  radiation  flux  at  the  thin  transition  boundary  layer  between  the  uniform 
plasma  core  and  the  ablative  surface;  a  model  of  the  transition  boundary  layer  to  obtain  the  boundary  conditions 
connecting  the  plasma  core  parameters  with  the  parameters  at  the  ablative  surface;  and  the  resistor-inductor- 
capacitor  circuit.  Thus,  the  model  allows  the  self-consistent  calculation  of  plasma  parameters  and  distribution  of 
wall  temperature  vs.  time.  The  radiation  wall  absorption  coefficient  is  assumed  to  be  constant  and  is  an  input 
parameter  of  the  model.  The  last  assumption  is  a  stretch  because  the  extinction  coefficient  of  polyethylene  is 
considerably  dependent  on  wavelength  and  temperature  which  may  change  significantly  in  the  course  of  discharge. 
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Since  the  model  is  not  capable  of  calculating  the  radiation  spectrum  at  the  ablative  surface  anyway,  this  is  the  extend 
of  what  model  is  capable  of. 

Our  model  shows  that  small  extinction  coefficients  of  the  wall  material  leads  to  large  heat  losses  from  the 
capillary  discharges  (the  heat  is  absorbed  by  the  bulk  of  the  capillary  wall  or  just  escapes  the  capillary),  to  the  spike 
in  the  plasma  temperature.  If  the  extinction  coefficient  is  too  small,  the  discharge  may  extinguish  because  the 
temperature  of  ablative  surface  does  not  increase  fast  enough  to  compensate  the  plasma  exhaust  from  the  open  end 
of  capillary.  Fig.  1.  Thus,  for  thruster  applications  the  capillary  wall  material  has  to  have  the  radiation  absorption 
coefficient  as  large  as  possible  to  provide  high  kinetic-to-electric-power  ratio. 


R 


AAA 


R, 


plasma 


Fig.  2  RLC  circuit  for  modeling  capillary  discharge. 


Slab  Capillary,  La  =  8  cm,  Da  =  4  mm 


Time  [s] 

Fig.  3.  LTE-factor:  blue  -  ij  =  0.5mm,  green  -  77  =  0.1mm,  red  -  77  =  0.02mm 
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Slab  Capillary,  La  =  8  cm,  Da  =  4  mm 


Time  [s] 

Fig.  4.  j}\  blue  -  i]=  0.5mm,  green  -  i]=  0.1mm,  red  -  77  =  0.02mm 


Slab  Capillary,  La  =  8  cm,  Da  =  4  mm 


Time  [s] 


Fig.  5.  Temperature  of  the  ablative  surface:  blue  -  77  =  0.5mm, 
green  -  77  =  0.1mm,  red  -  77  =  0.02mm 


Slab  Capillary,  La  =  8  cm,  Da  =  4  mm 


Time  [s] 

Fig.  6.  Polyethylene  number  density  in  the  plasma  core  region: 
blue  -  77  =  0.5mm,  green  -  77  =  0. 1mm,  red  -  77  =  0.02mm 


10 


Slab  Capillary,  La  =  8  cm,  Da  =  4  mm 


Time  [s] 

Fig.  7.  Plasma  pressure:  blue  -  77  =  0.5mm,  green  -  77  =  0.1mm,  red  -  77  =  0.02mm 


Slab  Capillary,  La  =  8  cm,  Da  =  4  mm 


Time  [s] 

Fig.  8.  Plasma  temperature:  blue  -  77  =  0.5mm,  green  -  77  =  0.1mm,  red  -  77  =  0.02mm 


Slab  Capillary,  La  =  8  cm,  Da  =  4  mm 


Time  [s] 

Fig-9.  Current:  blue  -  77  =  0.5mm,  green  -  77  =  0.1mm,  red  -  77  =  0.02mm 
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Slab  Capillary,  La  =  8  cm,  Da  =  4  mm 


Time  [s] 

Fig.  10.  Frad  and  0:  blue  and  light  blue  -  77  =  0.5mm,  green  and  light 
green  -  77  =  0. 1mm,  red  and  pink  -  77  =  0.02mm. 


Slab  Capillary,  La  =  8  cm,  Da  =  4  mm 


Fig.  11.  Wall  temperature  distribution:  77  =  0.02  pm,  blue  -  0  sec,  red  -  1 1  ps, 
green  -  1 1 1  ps,  and  brown  -  361  ps.  . 


Slab  Capillary,  La  =  8  cm,  Da  =  4  mm 


0.0  0.1  0.2  0.3  0.4 


x  [mm] 

Fig.  12.  Wall  temperature  distribution:  77  =  0.1  pm,  blue  -  0  sec,  red  -  1 1  ps, 
green  -  1 1 1  ps,  and  brown  -  361  ps.  . 
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Slab  Capillary,  La  =  8  cm,  Da  =  4  mm 


Fig.  13.  Wall  temperature  distribution:  ij  =  0.5  pm,  blue  -  0  sec,  red  -  1 1  ps,  green  -  1 1 1  ps,  brown  -  361  ps. 
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